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We have developed a semiclassical theory of short peri- 
odic orbits to obtain all quantum information of a bounded 
chaotic Hamiltonian system. If Ti is the period of the short- 
est periodic orbit, T2 the period of the next one and so on, 
the number Np,o of periodic orbits required in the calculation 
is such that Ti 4- ■ ■ ■ + TiVp „ ~ Th, with Th the Heisen- 
berg time. As a result Np.o — KTh / ^^ihTn), where h is the 
topological entropy. For methods related to the trace formula 



exp{hTH)/ihTH) 



PACS numbers: 05.45. +b, 03.65.Sq, 03.20.+i 
The semiclassical evaluation of energy spectra in clas- 
sically chaotic Hamiltonian systems was started in 1971 
with the Gutzwiller's trace formula From then on 
a big effort was dedicated to extend the formalism to 
wave functions j2| and to use resummation techniques to 
improve convergence properties of trace formulas 
However, a common drawback in all those approaches is 
the requirement of an enormous number of periodic or- 
bits (P.Os), restricting explicit calculations to very spe- 
cial systems; e.g., systems where classical mechanics is 
handled by a symbolic dynamics. Moreover, these ap- 
proaches miss a simple understanding of wave mechanics 
in terms of classical objects. 

Recently it was shown that chaotic eigenfunctions can 
be described in terms of a small number of localized struc- 
tures living on short periodic orbits Having this in 
mind, we present a simple semiclassical formalism to ob- 
tain quantum mechanics from a very small number of 
short P.Os. In order to do so, we will construct reso- 
nances (or scar functions) of short P.Os explicitly, 
evaluating the interaction between them. 

Let 7 be an unstable P.O isolated on each energy 
surface. We introduce a curvilinear coordinate system 
choosing the x axis along the trajectory and the y axis 
perpendicular to it at x (to simplify the exposition we 
take only one transversal direction). On the P.O y — 0. 
Classical mechanics in a neighborhood of the orbit is 
governed by a transversal symplectic matrix Al{x) of 
elements mij(x) = 1,2), which describes the lin- 
earized motion on the energy surface. Then, a point with 
transversal coordinates ( y, Py) at a; = evolves according 
to the following rule: 



Py{x) 



mii{x) mi2{x) 
m2i{x) m22{x) 
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We have selected the origin a; = such that mn (L) = 
(L) , with L the length of 7 . There are at least 2z^ 
points on the orbit satisfying that condition, with 1/ the 



Maslov index. With this choice the monodromy matrix 
M{L) acquires the form. 



M{L) = {-ly 



cosh(AL) sinh(Ai)/ tan((p) 
sinh(AL) tan((/3) cosh(AL) 



A is the Lyapunov exponent in units of [length^^] and 
tan((y9) (7^ 0) in units of [momentum /length] defines the 
slope of the unstable manifold in the plane y — Py (the 
slope of the stable manifold being — ta.n{(p)). 

In general it is impossible to compare vectors living in 
the plane y — Py because the axes have different units. 
However, when the directions are symmetrical with re- 
spect to the axes it is only necessary to compare one 
component. Then, we change to new axes ^„ and ^s, on 
the unstable and stable manifolds respectively, such that 
their projections on each axis are equal in absolute value. 
The symplectic matrix B transforming coordinates from 
the new axes into the old ones is 



B - (Cu 6) = il/V2) 



l/a ~s/a 



s a 



a 



with a = tan((p)| and s ~ sign{ip). Observe that 
B-^M{L)B = {-ly cxp{XLD) , with D a diagonal ma- 
trix of elements dn = 1 and ^22 = ^1- 

Now, we decompose M{x) into a periodic matrix F{x) 
describing the evolution of the manifolds, and a matrix 
(depending in a simple way of x) describing the exponen- 
tial contraction-dilation along the manifolds, 

M{x) = F{x) exp{xXK) = F{x) B exp(a;Ai:i) B-\ (1) 

with K = BDB-^ [^n /fc22 = 0, A:2i l/fci2 = tan((^)]. 
Equation (|l|) defines F{x) in terms of Al{x) and we can 
see that F{L) = (-1)"!. Floquet's theorem g affirms 
that the decomposition given in (^ can also be obtained 
for systems with many transversal directions. 

Now, it is possible to construct a family of resonances 
associated to 7. We adapt well known semiclassical 
techniques to obtain eigenfunctions concentrated in the 
neighborhood of stable P.Os In the x direction we 
consider the typical solution of a one dimensional motion, 
and in the transversal one we use a wave packet which 
evolves according to F{x), 



^yix,y) 



exp{i [S{x) + y2 r{x)/2]/n - i (l){x)/2} 



(2) 



where S{x) = mi dx is the action, T the period of 
7 and J the unit area in the plane y — py. T{x) = 



1 



P{x)/Q{x), with Q{x) [P{x)] the y (py) component of 
the complex vector 



(3) 



U^) +t = F{x){^^ +1 6) = M{x)B (^^J, 

Equation (||) shows that it is not necessary to evaluate 
explicitly F{x). The area preserving property of F{x) 
guarantees the following normalization condition: 

Q*{x) P{x) - Q{x) P*{x) = 2i A = 2iJ. (4) 

Then, Im[r(a;)] = J/\Q{x)\^ > 0. Accordingly V7 is 
concentrated around 7 and well-behaved in all regions 
except in the neighborhood of a turning point (i = 0) 
[ pT[ . Finally, the complex number Q{x) sweeps an angle 
(j^ix) while evolving from to x, and v = 0(i)/7r (this is 
the definition of the Maslov index). 

ipj is a continuous function at a; = L if the accumu- 
lated phase around the orbit is an integral multiple of 27r. 
This condition determines the admitted energies Ej of 
the P.O and corresponds to the Bohr-Sommerfeld quan- 
tization rule: S{L)/h — i't:/2 — 2mr, where n = 0, 1, . . . 
is the number of excitations along 7. 

We stress that the semiclassical construction of eigen- 
functions in the neighborhood of stable orbits is similar to 
. The initial complex vector of equation (||) is replaced 
by the eigenvector of the monodromy matrix (a complex 
vector in this case) satisfying (||). And of course, the evo- 
lution of the vector is given by the transversal symplectic 
matrix without modifications. Eigenvalues have an error 
0{h) and eigenfunctions an error O(v^). Moreover, it is 
possible to improve the accuracy by including transversal 
excitations Q. 

In our case, there is an essential error because the evo- 
lution in is given by a modified transversal matrix. 
In order to eliminate that error, we will first evaluate 
the action of the semiclassical evolution operator for in- 
finitesimal times over the resonance in the form 



H = ih limst^oiUiSt) - l)/6t. 



(5) 



The classical transversal evolution from x to x + 6x 
is given by M^{Sx) = M{x + 6x) M{xy^. So, 
Mx{5x) F{x) = F{x + 5x) exp{Sx XK). Observing that 
K^^—^y^ and K^g — ~^g^ it results to first order in dx: 



M^iSx) ^„(x) 



'{l + 5x X) £,u{x + 5x), 
(1 - 6x X) £,s{x + 5x). 



(6) 



The above expressions actually show clearly the approxi- 
mation involved in the construction. We have forced the 
vector {£,s) to evolve without dilation (contraction) 
while the right evolution dilates (contracts) the vector 
with a rate specified by A. Then, we set 5x ~ x 6t and 
use (^ to evaluate the action of U{6t) on -0^. Finally, 
after some calculations we obtain from iS], 



Hil)^{x,y) = g^{x,y) ijj^ix^y), (7) 
with g^{x, y) = E^+ihxX [y^J/h - \Q{x)\'^ /2)/Q{xf . 

Then, the application of the semiclassical Hamiltonian 
operator to the resonance gives the term E-ytp^ as ex- 
pected, plus a resonance of 7 with two excitations in 
the transversal direction. In fact, the two excitations 
are also expected because the right evolution produces a 
quadrupole-like deformation of the wave packet. Equa- 
tion (^ is a extremely powerful tool: it is the key to 
evaluate matrix elements between P.Os. As the oper- 
ator H is not exactly Hermitian for finite values of h, 
we define a symmetrized interaction between two P.Os 
7 and S as follows (in Dirac's notation): 



{6\H\j) ^ {{S\H^) 4 
{S\H^\j) ^ {H5\H^). 



{l\H5Y)/2, 



(8) 



By using (^, we obtain explicitly the following diagonal 
matrix elements in the semiclassical limit ?i — > 0, 



ii) 
Hi) 



(7l7> 



{l\H\l)/{l\l) ^E _ 
{l\H^\l)l{l\l)-E'^^{hXfx^/2^ 



7' 



with ±2 = S{L)/mT the time average of on the orbit. 
Expression Hi) shows that the width a-^ of the resonance 
is asymptoticaly proportional to A. Moreover, pE = 



o{h- 



|12(] shows that a unique orbit cannot support a 



stationary state in the semiclassical limit (to support an 
eigenfunction the width of a resonance needs to satisfy 
cr < 1). Of course, this result is well known Q. 

We should say some words about symmetry. If the sys- 
tem is time reversal it is easy to show that 'ip--y{x,y) = 
Tp^{x,y)* and Hip--f{x,y) — {Hip^{x,y))* , where —7 
is the time reversal partner of 7. If the system also in- 
cludes a spatial symmetry G, it results that ipcyix, y) = 
G-ip-yix^y) and H'il)G-y{x,y) — GHipj{x,y) . So, to ob- 
tain real eigenfunctions inside a defined symmetry repre- 
sentation, we construct real resonances inside the same 
representation by using group theory jl^ ]. 

For low or medium energies we can evaluate matrix 
elements directly on the domain; however, for high ener- 
gies or to obtain explicit expressions in terms of clas- 
sical quantities as fi goes to zero, it is preferable to 
work on a surface of section jl^. Let C be a differ- 
entiable curve with the coordinate q along it, and ry 
perpendicular to C at g (77 = on the curve). Sup- 
pose 7 crosses the section at (j = 1, . . . , m) with an- 
gles 9j] the corresponding positions on 7 being Xj (see 
Fig. 0(a)). In a neighborhood of radio 0{\/fi) around 
the intersection point j, the coordinates are related to 
order ?i by x — Xj = shi{9j) {q — qj) — cos{9j) rj and 
y ~ cos{9j) [q — qj) + sin(0j) 77 . Then, the restriction of 
"07 to up to 0{\/h)^ is a sum of wave packets in one 
dimension with tangential momentum pj = m Xj sin(0j). 
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(9) 



with i-)^^{q) ^ exp[i pj{q-qj)/h] ip^[xj,cos{ej){q-qj)]. 

If the system is bounded by a hard wall and we take 
this wall like the Poincare surface of section, -07 is null 
there up to order In the neighborhood of a 

bouncing point, ijjy consists of two terms associated to 
the incoming and outgoing trajectory such that the com- 
bination satisfies boundary conditions |^ (see Fig. 0(b)). 
Then, working with its normal derivative we define 



ifi di)^ 
2m drj 



{x,y)\^ =^ .i:jCos(6ij)?/'7,('?)- (10) 



Moreover, H(p^{q) = Hip^{x,y)\^ for (||) and Hip^{q) = 

-{ih/2m)dH^^/dr]{x,y)\^ for (|l^ are obtained to the 
leading order from (0) or (HG) taking into account that 



Hi^l, il) - 9'y[^],C0s{ej){q - qj)] ijj^^ (q). 



(11) 



From now on (p^ (q) (Eq. d) or (0)) is the object repre- 
senting the resonance on the section. In order to describe 
an effective Hilbert space we define a norm on 



(7l7>< 



(p.y{q)* (pj{q) f{q) dq, 



(12) 



such that i), ii) and Hi) are satisfied. Then, evaluating 
the leading term of the integral in (^2|), it results in a 
classical criterium for specifying f{q) 



(13) 



The signs (— ) and (-I-) correspond to (||) and ( p^ ) respec- 
tively. We select a smooth real function f{q), oscillating 
so slowly as possible, which satisfies ( p^ for all short 
P.Os required in the calculation. We notice that the 
existence of f{q) is not guaranteed for all sections. A 
hint to choose a section could be to make sure that the 
classical motion between consecutive points (the map) 
is simple. For example, in billiards the motion between 
bounces with the boundary is simple. On the other hand, 
eigenf unctions of billiards are reduced to the boundary in 
terms of their normal derivatives, and the metric on the 
boundary is defined by setting f{q) — 2 (r.n)(q)/i^ WA -, 
with n the unit outgoing normal to the boundary and r 
the position vector. Then, condition ( [l3| ) (with sign (-1-)) 
reduces to 

Em 
.^^ 2 (r.n)(,^)C0s(6lj), 

and this nice identity is valid for any P.O in any billiard; 
the demostration being a simple geometrical problem. 

Now we can evaluate matrix elements over the section 
by using equations (|^) or (pO|), (11) and ( p^ ) [S crosses 
the section at qk (A; = 1, . . . , m') with angles 9k], 



Ajk^l for 
Defining zi 



3-k=l 



A,k H^sM*0^iM) fi^)dq- (14) 



and Ajk — XjXk cos{9j) cos{6k) for (|l 
-icos^(6'i)r(a;i)/2, q = qi +ipi/2zi and 



Cjk = [zjCj + zlcD/lzj + z^), the Gaussian integrals in 
( |l4|) are given by 

D]k fjcjk) exp[-Bjk/h + i {aj - ak)] 
^T^ Ts Xj Xk \Q{xj) Q{xk)\ (zj + zl)/J 

Bjk = z.zlicj ~ clflizj + zl) + p]/Az, + pI/AzI and 
ai = S{xi)/h - (t){xi)/2. Moreover, Djk = 1 for O = i 
and Djk = g^[xj ,cos{9j)y/ {c^k - qjY + Ti/2{zj + zl)] for 
O = H . Then, Eq. ( p^ gives matrix elements in terms 
of classical quantities evaluated at the intersection points 
of the orbits with the surface of section. 

Finally, in order to obtain the eigenenergies and eigen- 
fuctions of a bounded chaotic Hamiltonian system in a 
given energy range, we proceed as follows. It is con- 
structed the family of resonances of the shortest periodic 
orbit 7i, living in the required energy range. The density 
of resonances associated to 71 is pi ~ Ti/2'!Th, with Ti 
the period of 71. Later, we do the same with the next 
shortest orbit 72, and so on (using only primitive orbits). 
The process stop when the whole density of resonances 
equals the mean energy density pg. 



k—1 ^ — 



(15) 



Equation ( [l5|) is actually impressive, it shows that the 
number of P.Os Np,o required in the calculation is very 
little and increases at most linearly with the Heisenberg 
time Th- More precisely Np,o — HTh / ^T^ihTn), where 
h is the topological entropy. For methods related to the 
trace formula Np,o — exp{hTH) / {hTn) ■ 

Another interesting quantity is the number of reso- 
nances Nres contributing to one eigenfunction. This 
number is proportional to the mean dispersion (sec Hi)) 
and to pe 



Nr 



3.6hpE v/< A2 ±2 >, 



(16) 



with < > the average over P.Os (using the factor 3.6 in 
(|l|), the 99% of an eigenfuntion is recovered). Then, we 
select N resonances (N > Nres), consecutive in energy, 
and call them Fi, . . . , Fat. Later, by solving the following 
generalized eigenvalue problem 



N 



^[{Tk\H\r,)-E{rk\r,))^,^o, \f k, (17) 



the eigenenergies E and eigenvectors ^ in the basis of 
resonances are given. 

The main idea for the selection of resonances is to 
obtain a quasi-orthogonal basis of highly localized (in 
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energy) functions. The best way of satisfying quasi- 
orthogonality is to use short periodic orbits. Now, for or- 
bits with comparable periods we select the one with mini- 
mum energy dispersion (see in)). This analysis works for 
hard chaos systems (where all P.Os are unstable and iso- 
lated). However, for systems with a fraction of regular 
motion, we need to include the same fraction of regular 
functions in the basis. And, for systems with a continu- 
ous family of neutral P.Os (e.g. the bouncing ball family 
in the stadium billiard), it is required the corresponding 
fraction of phase space localized functions. 

In conclusion we need to construct adequate functions 
in each classically different region of phase space, the 
number of them satisfying the required mean density (ob- 
tained semiclassically) in each region. In a chaotic region, 
functions (we call them resonances) are constructed with 
the shortest P.Os, and the number of them used to fill 
the Hilbert space is very little as implied in Eq. (p^. In 
order to obtain the interaction between two given short 
P.Os it is possible to follow different strategies. Think- 
ing at classical level it is necessary to use at least a P.O 
living in the neighborhood of the previous ones. Then, 
to obtain all matrix elements, the period of the orbits 
required in the full calculation would be of the order of 
the Heisenberg time, and no advantage is reached with 
respect to other approaches. However, in this article we 
showed that thinking at quantum level, the interaction 
between short P.Os can be evaluated simply in terms of 
transversal excitations. In this way, all the information 
required in the calculation is contained in short P.Os. 

We finish with some remarks, i) This theory has been 
applied successfully to the desymmetrized stadium bil- 
liard ii) In this theory, eigenvalues have an error 
0{h) and eigenfunctions an error 0{Vh). in) We be- 
lieve that the inclusion of transversal excitations in the 
construction of resonances is the way of extending this 
theory to higher orders in h. iv) The basis of resonances 
is particularly useful for parametric dependent systems. 
In fact, this is the diabatic basis, so difficult to find in 
chaotic systems 1^. v) The evaluation of func- 
tion of n is a simple an efficient semiclassical measure 
of scars . vi) The phenomenon of chaos-assisted tun- 
nelling can be studied without external parameters [po| . 
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